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Abstract 

We show that through an interdot off-site electron correlation in a double 

quantum-dot (DQD) device, Kondo resonances emerge in the local density of 

states without the electron spin-degree of freedom. We identify the physical 

mechanism behind this phenomenon: rather than forming a spin singlet in the 

device as required in the conventional Kondo physics, we found that exchange 

of electron position between the two quantum dots, together with the off-site 

Coulomb interaction, are sufficient to induce the Kondo resonance. Due to 

peculiar co-tunneling events, the Kondo resonance in one dot may be pinned 

at the chemical potential of the other one. 
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The recent discovery of Kondo effect in semiconductor quantum dots (QD) |],[| has 
generated tremendous interests both theoretically and experimentally When a QD is con- 
tacted by two metallic leads, Kondo effect arises as a result of coherent superposition of 
co-tunneling processes |2|U such that, at low temperatures, a very narrow Kondo resonance 
emerges in the density of states of the QD at the free electron chemical potential /i, and a 
spin singlet is formed in the whole system of the QD plus the leads. Kondo effect has so 
far been investigated in many systems and its understanding provided important insights to 
strongly interacting electrons in nanostructures. In most previous theoretical investigations, 
the physical systems under consideration all possess three common and fundamental ingre- 
dients: [U the existence of co-tunneling (i.e. exchange) events, the role of the electron spin 
degrees of freedom, and finally the on-site Coulomb interaction U Q . Because the QD Kondo 
effect is a result of many-body correlation, the existence of a great deal of co-tunneling 
events and the Coulomb interaction are crucial. However, must it be the on-site interac- 
tion? Are spin degrees of freedom essential? Although no exception has been found so far, 
these questions are intriguing because their understanding will provide a clear picture to 
the deciding physical factors of QD Kondo phenomenon. Our investigation found that QD 
Kondo effect does occur even without considering the spin degrees of freedom and without 
the on-site Coulomb energy — provided that there are spatial degrees of freedom and there 
is an off-site interaction. It is the purpose of this letter to report our investigations on these 
issues. 

In particular, we investigate a two-probe double QD (DQD) device in which the spin 
degeneracy is neglected. The two QDs provide a spatial degree of freedom for co-tunneling 
processes (see below) and there is an off-site e-e interaction between the dots. It should 
be mentioned that conventional Kondo effect, i.e. the Kondo effect caused by the three 
fundamental ingredients mentioned in the last paragraph, in DQD systems have been the 
subject of several recent studies 0-0. These investigations considered effects of the intradot 
on-site interaction U and/or interdot magnetic spin exchange. The latter is not an inde- 
pendent effect as it is determined by the (intradot) interaction U and the interdot tunnel 
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coupling [0. The present work, however, focuses on Kondo effect induced by the interdot 
off-site interaction U without the spin degrees of freedom. Our main findings are: (i) At 
low temperatures there are Kondo resonances in the local density of states (LDOS) of the 
DQD which occur without the entire system forming a spin singlet; (ii) If the interdot tun- 
nel coupling is zero, the Kondo resonance in one QD is pinned at the chemical potential of 
the lead that is attached to the other dot; (iii) At a non-zero interdot tunnel coupling, this 
Kondo resonance is split so that sharp features including peaks and abruptly drops emerge 
in the conductance of the device. These results are significant because they show that the 
Kondo resonance can occur by an off-site interaction even without the electron spin-degrees 
of freedom. 

Our system of DQD coupled with two leads is described by the following Hamilto- 
nian, H = H D q D + H T + H a , where H a = e ak< 1 ak a ak describe the noninteraction 

a(a=L,R) k 

left/right lead; Ht = Y^k^a^ak^a + H.c.) denotes the tunneling between the left (right) 
QD and the left (right) lead. Hdqd = J2 a e a<$ a d a + Ud L dLd R dn + tc{d\dR + H.c.) models 
the coupled double QD, in which each QD has only a single level without spin index. The 
purpose of neglecting the spin index is to emphasize that fact that for this system, Kondo 
effect occurs without the need of spin-degrees of freedom. Realistically, if a high magnetic 
field is applied to the QDs, leading the spin separation being larger than the QD energy level 
spacing ||, the opposite spin state essentially does not affect tunneling at low bias because 
its contribution is actually weaker than that of the first excited QD state. In this situation 
the spin-degrees of freedom can be neglected. Finally, when bias voltage V is not very high, 
only one eigenstate on each QD is active, therefore the on-site interaction U of each QD 
can be absorbed into its single particle level e a (a = L,R), || and only the interdot off-site 
interaction U will be considered. 

Our theoretical analysis is from the Keldysh nonequilibrium Green's functions. The 
current from the left lead flowing into the left QD can be expressed as (h = 1) [[K]]: I = 
—2elm J L{fL(c)G r LL (e) + |G^ L (e)}, where f a (e) is the Fermi distribution of the a lead, 
and T a = 2nJ2k |^a| 2( K e — e fc)- The Green's function G r £ Ij (e) are the Fourier transformation 
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of G r it(t), and G r LL {t) = -i6(t) < {d L (t), 4(0)} >, G< L (t) = % < d{(0)d L (t) >. Using the 
standard equation of motion technique and taking the familiar decoupling approximation 



I!]] , we have solved G r (e) in the infinite U limit (U — ► oo) to be: 
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where T, d = + Y, R . S^(e) = J2k \ta\ 2 / (e — e«fc + is the lowest-order self-energy for a 
noninteraction system; the higher-order self-energy , are: (e) = Z)fc[(e — e a fc)(e^= 

Ae-e afc )-2|t c | 2 ] J B Qfc (e), E£(e) = £ fc 2|t c | 2 J B afc (e), and E£(e) = £ fc (e T Ae - e ak )t c B ak (e), 
withB ak (e) = \t a \ 2 f a (e ak )/(e-e ak )/[(e-Ae-e ak )(e+Ae-e ak )-4:\tc\ 2 ] and Ae = e L -e R . The 
quantity n a =< d) a d a > in Eq.(l) is the intradot electron occupation number, < d) L d R > 
and < d R di, > are the interdot correlation due to the interdot tunnel coupling and the 
interaction. They must be calculated self-consistently: < d^dp >= —i J ^G^ a {e). It should 
be mentioned that although the equation of motion method only gives a correct qualitative 
physics |[L2|| , it is sufficient for the purpose here. 

Next we solve the Keldysh Green's function G < (e). For interacting systems, G < (e) cannot 
be obtained from equation of motion without introducing additional assumptions flT2"| , |T3 



In previous work, various approximations were developed including the use of an ansatz for 
interacting lesser and greater self-energy ||13|| : the noncrossing approximation f[2" |, and so 
on. Any of them could be sufficient for our purposes. However, we note that since only 
J deG < (e) is actually needed for calculating both the electric current and for iterating the 
self-consistent equation, one does not need to solve G < (e) itself. We found that quantity 
/ dtG K (e) can actually be solved exactly, allowing us to bypass any approximation involved 
in computing G < (e). We now present this technical advance which is quite general and can 
be applied to many other systems. As stated above, J deG^ a (e) ~< d^dp >, we therefore 
proceed to first derive the operator $ a (f)dp(f) (a, (3 = L,R) over time t, and then take 
average < 4;[dl [ (t)dp(t)] >= 0. The last equality is true because we consider steady state 



transport. On the other hand, the time derivative can be explicitly calculated using the 
Heisenberg equation, and we obtain 

(e a - e Q - i — ^ Tr ) < da > +t c < S a d a > -t c < d^ >= J ^T a f a G r aa - J ^T a f a G% c 

(2) 

-t c < did,, > +t c < d\d a > -iT & < d)^ >= f p-TMG^ - G%^\ (3) 

J ATI 

where a = R for a = L, or a = L for a = R. From Eqs.(^|,H) and having already 
solved G r (e), < d^dp > (i.e. J deG^ a (e)) can be exactly solved and current can now be 
obtained without further difficulty In our numerical evaluations, we assume square bands 
of width 2W and symmetric coupling barriers so that r^(e) = r#(e) = T9{W — |e|), with 
W = 1000 3> max(kBT, eV, T). Here V = Hl — is the bias voltage. We also fix Y = 1 as 
energy unit and let hl = —j^R- 



First, we investigate LDOS of the DQD device. [|J] Because the LDOS of the right dot is 
similar to that of the left one, we only discuss the latter. Consider equilibrium situation (bias 
V = 0) and the case where the two dots have equal level position (Ae = — e# = 0), the 
LDOS is shown as the right inset of Fig. Id in which various curves correspond to different 
values of the interdot tunnel coupling tc- When there is no tunnel coupling (Tc = 0, thin 
solid line), besides a broad main peak at e ~ —1.5 (not shown) which corresponds to the 
renormalized left dot level e^, a sharp peak emerges at energy e = \i (now /i = = ixr). 
This peak has the typical Kondo characteristic: it only exists at low temperature and is 
pinned at e = ji. Turning on the tunnel coupling (t c ^ 0) between the two dots, the Kondo 



peak is split into two smaller peaks at e = /i ± 2tc and whose weight are \ |15| . Setting 



the off-site Coulomb energy U = 0, the Kondo peaks disappear. Clearly, it is the interdot 
off-site Coulomb interaction U that induces the Kondo peak. 

The origin of the Kondo peaks at tc = is from an interesting co-tunneling process 
which is shown in the left inset of Fig. Id. When e^, €r < Hl, Hr < £l + U,€r + U, the DQD 
has only one electron (e.g. it occupies the right dot) and the first order tunneling is blocked 



by the off-site interaction. However, due to Heisenberg uncertainty principle, higher-order 
co-tunneling processes consisting of two virtual tunneling events can still take place. In such 
a process, the electron in the right dot tunnels to the Fermi level of the right lead, followed by 
an electron in the left lead at the Fermi level tunneling into the left dot, on a very short time 
scale ~ %/{hr — This co-tunneling process, in effect, moved the original electron from 
the right dot to the left one. At low temperatures, a frequent exchange of electrons between 
the two dots occur through this process even though tc = 0, and a coherent superposition 
of all possible co-tunneling processes of this type produces the Kondo resonance discussed 
in the last paragraph in which a very narrow peak emerges at e = \i in the LDOS curve 
of left (right) dot. We emphasize that this mechanism of inducing a Kondo resonance is 
qualitatively different from the co-tunneling process that induces the conventional Kondo 
effect by the on-site interaction: in that case the local spin is flipped after the co-tunneling 
process. 

Next, we discuss the non-equilibrium situation (V ^ 0) keeping Ae = 0. Fig. Id shows 
the left dot's LDOS. At tc = 0, a narrow Kondo peak emerges (thin solid line). However, 
its position is at e = fiR, not at In other words, the Kondo peak of one dot occurs at 
the chemical potential of the lead of the other dot! This is a very surprising result indeed. 
It is also qualitivally different from the Kondo resonance induced by intradot interaction U Q 
(not the interdot U) where the peak is pinned to the chemical potential of the lead attached 
to the same dot f|. When tc ^ 0, states and €r hyberdize into two "molecular" states 
e ± = £i±£a ± A^ expanding to the entire device at Ae = fL6|, where AE = ^ Ae 2 + At 2 c . 
The corresponding left dot's LDOS are plotted in Fig. Id. The main features are: (i) The 
Kondo peak at e = /j,r when tc = is now split into three peaks, at e = fiR and //r ± AE 
with weights | and | respectively, (ii) A dip emerges at e = \i L and two more peaks are at 
fii ± AE. Most importantly, all these (five) peaks and the dip have Kondo characteristics: 
(1) they are pinned at (jll/r and hl/r± AE; (2) they disappear at high temperature; on the 
other hand, at lower temperature they become higher/lower for the peak/dip (See Fig.2); 
(3) they all disappear when U = 0. 
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In fact, the five Kondo peaks and the dip are results of many-body correlation due to the 
off-site interaction U. They originate from the co-tunneling processes between the molecular 
state e ± and the two leads which are shown in Fig.la-c. For example, Fig. la describes the 
following co-tunneling process. An electron initially occupies e + and tunnels to the Fermi 
level of the right lead. At almost the same moment another electron in the left lead, having 
energy hr — AE, tunnels into the other molecular state e~ . At low temperatures, a coherent 
superposition of many co-tunneling events of this type gives rise to the Kondo resonance at 
fi R — AE in the LDOS. Exactly the same way, the co-tunneling events shown in Fig.lc induce 
a Kondo peak at hr + AE. Fig. lb shows a different type co-tunneling process. Here an 
electron occupies level e ± initially. It then tunnels to the Fermi level of the right lead, closely 
followed by another electron in the left lead with energy /i R which tunnels into same level e ± . 
This co-tunneling process does not induce electron exchange between the molecular states 
e^, therefore it cannot induce a Kondo peak at the DOS of the molecular states. However, 
molecular states are linear combinations of the original dot states 6l and €r, so in this 
co-tunneling process, electrons may still exchange between e L and e R , which gives rise to a 
Kondo peak at e = /ir in the left dot's LDOS and a dip at e = /ir in the right dot's LDOS. 

The peculiar Kondo resonances discussed so far are for equal dot states level positions, 
i.e. for Ae = €l — £r = 0. They however persist even when Ae ^ 0. As Ae increases from 
zero, the molecular states e ± gradually localize to one of the dots and they asymmetrically 
couple to the two leads. If Ae > 0, e + /~ are coupled stronger to the left/right lead, and 
weaker to the other lead (Fig.la-c). Therefore, the co-tunneling processes of Fig.(la,lb) will 
be weakened while that of Fig.lc enhanced. Correspondingly, in the left dot's LDOS, the 
Kondo peaks at /j,r + AE is increased, and others will be reduced (see inset of Fig. 2). In 
the limiting case of t c j Ae — > 0, only the Kondo peak at fi R + AE survives in the left dot's 
LDOS. It is not difficult to find that approximately the weights and positions for each left 
dot's Kondo peak are: |(1 — — H^Ja for Kondo peak at fi R — AE; for Kondo peak 
at fJ'R'i |(1 + x§) — for Kondo peak at hr + AE; j^^L for Kondo peak at hl — AE; 
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~ A £2 for Kondo dip at fiL] and for Kondo peak at hl + AE. 

The current I and differential conductance dl/dV are substantially inffuenced by the 
Kondo resonances. The left inset of Fig. 3 show dl/dV versus V for several different in- 
terdot coupling tc- Since we fixed Ae = 0, dl/dV must be symmetric while / must be 
antisymmetric across zero bias which is what we found. The main characteristics are: (1) 
A very narrow peak is exhibited at V = in dl/dV, as a result of the co-tunneling process 
of Fig. lb. Increasing tc, the height of this peak first increases sharply and then reaches a 
maximum ~ O.Qe 2 /h (the unitary value is e 2 /h) at t c ~ 0.18r, followed by a slow decrease. 
(2) dl/dV abruptly drops at V = ±AE. If tc is large, dl/dV exhibits another peak at 

V = ±AE followed by the drop. Because at V = ±AE, the left lead's chemical potential 
Hi, aligns with the left dot's Kondo resonance at /ir ± AE, while /ir aligns with the right 
dot's Kondo resonance at /il T AE. (3) Increasing temperature, the dl/dV features at 

V = and V = ±AE gradually diminish, and the overall value of dl/dV greatly reduces 
(see the right inset of Fig. 3). At high temperature, e.g. T = 0.5, dl/dV becomes very small 
and very weakly dependent on bias. 

When Ae 7^ 0, dl/dV is no longer perfectly symmetric as shown in Fig.3. With increasing 
Ae, electron tunneling between the left and right dots becomes more difficult and molecular 
states e ± will localize to one of the dots. As a consequence, the features of dl/dV at 

V = and ±AE are reduced. Specifically, for Ae > the abrupt drop of dl/dV at 

V = —AE will gradually disappear, because the co-tunneling process of Fig. la becomes 
weaker. On the other hand, the abrupt drop at V — AE survives although it becomes 
weaker. Furthermore, as discussed above, Ae > enhances the co-tunneling process of 
Fig.lc, due to this enhancement a new peak in dl/dV arises at V — AE for large Ae. 

In summary, we investigated Kondo effect in a DQD device with an off-site electron 
correlation provided by the interdot interaction U . A number of co-tunneling processes have 
been identified between the molecular states and the leads, and electrons may be moved 
from one dot to the other due to these processes. Rather than flipping the electron spin, the 
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interdot exchange of electrons together with the off-site electron correlation and coherent 
superpositions of the co-tunneling processes, induce several Kondo resonances without the 
electron spin-degrees of freedom. We note that although the intradot interaction U Q should 
be larger than U, the Kondo temperature is not monotonic in these parameters. It is 
actually quite probable that the Kondo resonance determined by U is more prominent than 
that induced by U Q . Therefore, if spin, intradot and interdot interactions are all present, 
more complicated Kondo phenomenon should result. 
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FIGURES 

FIG. 1. (a)-(c) are schematic diagrams of co-tunneling processes between the two molecular 
states and the two leads, (d) and its right inset show the left dot's LDOS vs energy e at different 
tunnel coupling tc with e^ = €r = —2.5, T = 0.005. In (d), [1l = —^r = 0.2; in its right inset 
/j,L = —/~ir = 0. The left inset of (d) is a schematic diagram for the co-tunneling process at tc = 0. 

FIG. 2. Left dot's LDOS vs e at different temperature T. Inset: data for different Ae and 
£l/zr = —2.5 ± Ae/2. tc = 0.05 and other parameters are the same as those of Fig. (Id). 

FIG. 3. The conductance dl/dV vs the bias V at different Ae. el = —2, cr = — 2 — Ae, 
T = 0.001, and t c = 0.05. Left inset: dl/dV vs V at different t c with e L = e R = -2 and 
T = 0.001. Different curves correspond to tc = 0.025, 0.05, 0.1, and 0.2 from bottom to top. Right 
inset: dl/dV vs V at different temperature T where tc = 0.05 and e^ = €r = —2. 



12 




Fig.2 



